Fourier spectral methods for fractional-in-space reaction-diffusion equations

Fractional differential equations are becoming increasingly used as a powerful modelling approach for understanding the many aspects of nonlocality and spatial heterogeneity. However, the numerical approximation of these models is computationally demanding and imposes a number of computational constraints. In this paper, we introduce Fourier spectral methods as an attractive and easy-to-code alternative for the integration of fractional-in-space reactiondiffusion equations. The main advantages of the proposed schemes is that they yield a fully diagonal representation of the fractional operator, with increased accuracy and efficiency when compared to low-order counterparts, and a completely straightforward extension to two and three spatial dimensions. Our approach is show-cased by solving several problems of practical interest, including the fractional Allen–Cahn, FitzHugh–Nagumo and Gray–Scott models,together with an analysis of the properties of these systems in terms of the fractional power of the underlying Laplacian operator

Simulation of cell movement through evolving environment: a fictitious domain approach

A numerical method for simulating the movement of unicellular organisms which respond to chemical signals is presented. Cells are modelled as objects of finite size while the extracellular space is described by reaction-diffusion partial differential equations. This modular simulation allows the implementation of different models at the different scales encountered in cell biology and couples them in one single framework. The global computational cost is contained thanks to the use of the fictitious domain method for finite elements, allowing the efficient solve of partial differential equations in moving domains. Finally, a mixed formulation is adopted in order to better monitor the flux of chemicals, specifically at the interface between the cells and the extracellular domain

Evaluation of the Gulf Fisheye Bycatch Reduction Device in the Northern Gulf Inshore Shrimp Fishery

The performance of the Gulf fisheye bycatch reduction device (BRD) was evaluated on two vessels during inshore shrimp fishing operations in the northern Gulf of Mexico by comparing catch rates with control nets in twin-trawl configurations using typical inshore nets with 7.6-m headropes. The BRD produced substantial reductions in finfish bycatch with no shrimp (Penaeus spp.) loss in three of the four evaluations. Proper installation of the BRD in the net is critical in maximizing bycatch reduction and preserving the shrimp catch. For the inshore fishery in the northern Gulf of Mexico, the recommended distance to install the fisheye BRD in front of the bag tie on 7.6-m headrope shrimp trawls is 2.6 m. Because no shrimp loss attributed to the BRD was noted during all four evaluations in this study, the results suggest that the Gulf fisheye BRD could be used effectively year-round in the northern Gulf inshore fishery and not just when finfish are overly abundant on the shrimp fishing grounds

'Extremotaxis': Computing with a bacterial-inspired algorithm

We present a general-purpose optimization algorithm inspired by “run-and-tumble”, the biased random walk chemotactic swimming strategy used by the bacterium Escherichia coli to locate regions of high nutrient concentration The method uses particles (corresponding to bacteria) that swim through the variable space (corresponding to the attractant concentration profile). By constantly performing temporal comparisons, the particles drift towards the minimum or maximum of the function of interest. We illustrate the use of our method with four examples. We also present a discrete version of the algorithm. The new algorithm is expected to be useful in combinatorial optimization problems involving many variables, where the functional landscape is apparently stochastic and has local minima, but preserves some derivative structure at intermediate scales

Fractional diffusion models of cardiac electrical propagation: role of structural heterogeneity in dispersion of repolarization

Structural heterogeneity constitutes one of the main substrates influencing impulse propagation in living tissues. In cardiac muscle, improved understanding on its role is key to advancing our interpretation of cell-to-cell coupling, and how tissue structure modulates electrical propagation and arrhythmogenesis in the intact and diseased heart. We propose fractional diffusion models as a novel mathematical description of structurally heterogeneous excitable media, as a mean of representing the modulation of the total electric field by the secondary electrical sources associated with tissue inhomogeneities. Our results, validated against in-vivo human recordings and experimental data of different animal species, indicate that structural heterogeneity underlies many relevant characteristics of cardiac propagation, including the shortening of action potential duration along the activation pathway, and the progressive modulation by premature beats of spatial patterns of dispersion of repolarization. The proposed approach may also have important implications in other research fields involving excitable complex media

Runge-Kutta methods for third order weak approximation of SDEs with multidimensional additive noise

A new class of third order Runge-Kutta methods for stochastic differential equations with additive noise is introduced. In contrast to Platen's method, which to the knowledge of the author has been up to now the only known third order Runge-Kutta scheme for weak approximation, the new class of methods affords less random variable evaluations and is also applicable to SDEs with multidimensional noise. Order conditions up to order three are calculated and coefficients of a four stage third order method are given. This method has deterministic order four and minimized error constants, and needs in addition less function evaluations than the method of Platen. Applied to some examples, the new method is compared numerically with Platen's method and some well known second order methods and yields very promising results.Comment: Two further examples added, small correction

DBI Galileon and Late time acceleration of the universe

We consider 1+3 dimensional maximally symmetric Minkowski brane embedded in a 1+4 dimensional maximally symmetric Minkowski background. The resulting 1+3 dimensional effective field theory is of DBI (Dirac-Born-Infeld) Galileon type. We use this model to study the late time acceleration of the universe. We study the deviation of the model from the concordance \Lambda CDM behaviour. Finally we put constraints on the model parameters using various observational data.Comment: 16 pages, 7 eps figures, Latex Style, new references added, corrected missing reference

de Sitter Galileon

We generalize the Galileon symmetry and its relativistic extension to a de Sitter background. This is made possible by studying a probe-brane in a flat five-dimensional bulk using a de Sitter slicing. The generalized Lovelock invariants induced on the probe brane enjoy the induced Poincar\'e symmetry inherited from the bulk, while living on a de Sitter geometry. The non-relativistic limit of these invariants naturally maintain a generalized Galileon symmetry around de Sitter while being free of ghost-like pathologies. We comment briefly on the cosmology of these models and the extension to the AdS symmetry as well as generic FRW backgrounds

Discrete Razumikhin-type technique and stability of the Euler-Maruyama method to stochastic functional differential equations

A discrete stochastic Razumikhin-type theorem is established to investigate whether the Euler--Maruyama (EM) scheme can reproduce the moment exponential stability of exact solutions of stochastic functional differential equations (SFDEs). In addition, the Chebyshev inequality and the Borel-Cantelli lemma are applied to show the almost sure stability of the EM approximate solutions of SFDEs. To show our idea clearly, these results are used to discuss stability of numerical solutions of two classes of special SFDEs, including stochastic delay differential equations (SDDEs) with variable delay and stochastically perturbed equations